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ABSTRACT 



The POINT-AGAPE collaboration is currently searching for massive compact halo objects (MA- 
CHOs) towards the Andromeda galaxy (M31). The survey aims to exploit the high inclination of the 
M31 disk, which causes an asymmetry in the spatial distribution of M31 MACHOs. Here, we investigate 
the effects of halo velocity anisotropy and flattening on the asymmetry signal using simple halo models. 
For a spherically symmetric and isotropic halo, we find that the underlying pixel-lensing rate in far-disk 
M31 MACHOs is more than 5 times the rate of near-disk events. We find that the asymmetry is increased 
further by about 30% if the MACHOs occupy radial orbits rather than tangential orbits, but is substantially 
reduced if the MACHOs lie in a flattened halo. However, even for haloes with a minor-to-major axis ratio 
q = 0.3, the numbers of M31 MACHOs in the far-side outnumber those in the near-side by a factor of 
PsJ . ^2. There is also a distance asymmetry, in that the events on the far-side are typically further from the 

major axis. We show that, if this positional information is exploited in addition to number counts, then the 
number of candidate events required to confirm asymmetry for a range of flattened and anisotropic halo 
. models is achievable, even with significant contamination by variable stars and foreground microlens- 

\ ing events. For pixel-lensing surveys which probe a representative portion of the M31 disk, a sample of 

pi i' around 50 candidates is likely to be sufficient to detect asymmetry within spherical haloes, even if half the 

Q . sample is contaminated, or to detect asymmetry in haloes as flat as ^ = 0.3 provided less than a third of 

)—{ \ the sample comprises contaminants. We also argue that, provided its mass-to-light ratio is less than 100, 

the recently observed stellar stream around M3 1 is not problematic for the detection of asymmetry. 

Subject headings: Dark Matter - Galaxies: Individual (M3 1) - Gravitational Lensing 
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1. Introduction 

Continuing disagreement as to whether Massive Compact Halo Objects (MACHOs) have been detected by mi- 
crolensing experiments looking towards the Magellanic Clouds highlights the need for other microlensing targets (e.g., 
Kerins 2001). The Andromeda Galaxy (M31) presents an opportune target in this respect. The disk of M31 is highly 
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inclined (i ~ 77°), with the consequence that Hnes of sight to disk stars in the north-west or near side of M3 1 are shorter 
than those to the south-east or far side. Microlensing by a spheroidal dark halo will have a characteristic signature with 
an excess of events on the far side of the M31 disk (Crotts 1992; Baillon et al. 1993). This asymmetric signal is absent 
for variable stars or stellar microlenses in the disk of M31. A number of groups (e.g., Auriere et al. 2001; Riffeser 
et al. 2001; Calchi-Novati et al. 2002; Crotts et al. 2001) are now carrying out large-scale surveys of M31 to search 
for this near- far disk asymmetry. This is a mammoth task as the individual stars in M3 1 are not resolved, so that new 
techniques based on the super-pixel method (Ansari et al. 1997) or difference imaging (e.g., Crotts & Tomaney 1996) 
have been exploited to measure the flux changes on unresolved stars. Nonetheless, convincing candidate events are 
now being discovered, for example by the POINT- AGAPE collaboration (e.g., Auriere et al. 2001; Paulin-Henriksson 
et al. 2002, 2003). Therefore, this is a timely moment to consider what factors affect the near-far disk asymmetry and 
how many events are likely to be needed for a convincing detection. 

The aim of this paper is to estimate the size of candidate event samples needed to detect asymmetry for different 
halo models. In Sections 2 and 3 of the paper we show how the magnitude of the asymmetry signal is affected by the 
velocity anisotropy and the flattening of the M3 1 baryonic dark halo respectively. There are few ways known to us 
for measuring the properties of the orbits of dark objects in the halo of any galaxy or for ascertaining the flattening of 
the baryonic dark component of the halo. Hence, any clues gleaned from pixel lensing experiments will be invaluable. 
In Section 4, we present some simple, non-parametric statistical estimators of asymmetry and calculate how many 
candidates are needed to give convincing detections. 



2. The Effects of Velocity Anisotropy 
2.1. Anisotropic Models 

Here, we investigate the microlensing properties of haloes in which the velocity distribution of the MACHOs is 
anisotropic. We use models in which the halo density p is isothermal and the rotation curve is flat: 

poc-^, Vn = constant. (1) 

For these models, the velocity dispersions in the spherical polar coordinate system are given by (see e.g.. White 1981; 
Evans, Hafner & de Zeeuw 1997) 

al = al = il + a)aj=^-^, (2) 

where a > -1 is the anisotropy parameter. If a < 0, then the velocity distribution is referred to as 'radially anisotropic' ; 
if a > 0, then it is 'tangentially anisotropic' Whilst ct^ inequation (2) diverges as a — »-l, velocity dispersion ratios are 
rarely observed to be more extreme than 3:1, implying -8/9 < a < 8 . We assume a circular velocity of vq = 235 km s"' 
for M3 1 's halo (Emerson 1976) and we compute the cutoff radius to give a total halo mass of 1 x 10^^ M© (e.g., Evans 
& Wilkinson 2000). The sources are drawn from the M3 1 disk which is adequately modelled as a sheet inclined at 77° 
to the line of sight. The source velocity is assumed to be dominated by the disk rotation speed of 235 km s~^ 

For each halo model, we calculate a theoretical estimate of the pixel-lensing rate Fp. Unlike the classical (resolved 
star) microlensing rate, Fp depends additionally on the surface brightness of the M3 1 disk and the luminosity function 
of the M3 1 sources (c.f., Kerins et al. 2001). The calculations here are performed for a V-band luminosity function and 
surface brightness distribution. The surface brightness of M31 is tabulated in Walterbos & Kennicutt (1987). The M31 
disk luminosity function is assumed to be the same as for the Milky Way and we use the data of Wielen, JahreiB & 
Kruger (1983) to characterise the faint end (My > 5) and that of Bahcall & Soneira (1980) for the bright end. For each 
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source star of flux F at sky position (x,y) on the M3 1 disk, we compute the maximum impact parameter, uj, needed 
to ensure that the magnified source star noticeably enhances the local flux contribution from the background galaxy 
and sky background. Motivated by the POINT-AGAPE experiment (e.g., Paulin-Henriksson et al. 2002, 2003), we 
assume that detection is performed on a "super-pixel" array of pixels of size 2."1 x 2."1 which typically encloses 40% 
of the total flux from a point source. In addition to the background due to the M31 surface brightness we allow for a 
sky background of 19.5 mag arcsec"^. We assume an event is detectable if the flux change caused by microlensing is 
> 1 % of the background flux on the super-pixel. The pixel lensing rate is then 

TpCx,}') = (MT(x,>'))rc(x,>'), (3) 

where Fc is the classical microlensing rate (Paczynski 1986; Kiraga & Paczynski 1994) and {iij) is the maximum 
impact parameter averaged over the luminosity function </> of the source stars. Explicitly, we can write that 

This is an upper limit to the observed pixel lensing rate for any real experiment, as it does not take into account the 
effects of sampling, changing observing conditions or event identification algorithms. These effects may alter the 
spatial distribution of events, but they can be corrected via the calculated detection efficiency. In fact, the efficiency 
is largely controlled by the local surface brightness and so is approximately symmetric with respect to the major axis. 
The ratio of the number of far-disk to near-disk events therefore does not depend on the efficiency to lowest order 

Finally, the spatially-averaged pixel-lensing rate (Fp) is obtained by weighting the rate with the M3 1 disk surface 
brightness. The central portions of the M3 1 disk are omitted, partly because stellar lenses in the M3 1 bulge dominate 
here and partly because the halo model is singular at the centre. So, a central region of 5' radius is excised from the 
M31 disk before performing the spatial averaging separately for events above and below the M31 major axis. None 
of the current experiments is surveying the entire M3 1 disk, but their fields do span the large majority of the minor 
axis and we therefore expect any underlying asymmetry in their fields to be representative of the globally-averaged 
asymmetry computed here. 



2,2. Results 

Figure 1 shows the spatially-averaged theoretical pixel-lensing rate, (Fp), for lenses of mass M, normalised to 
the value Fq = 7.6 x 10"^ (M/M©)"'/^ events per star per year. When a < 0, the velocity distribution is radially 
anisotropic and the rate (Fp) oc vo/Vl + a. When a > 0, the velocity distribution is tangentially anisotropic and the 
rate (Fp) cx vo, so the total rate in radially anisotropic models is somewhat higher. However, we are primarily interested 
in the differences between the near and far disk. Such asymmetries may manifest themselves in the numbers, locations 
and time-scales of the events. Accordingly, Figure 1 also shows the far-to-near disk ratios for the pixel-lensing rates 
(A = (rp)f/ (rp)n), the mean Einstein crossing times {{tn)tl (^E)n) and the ratio of mean projected distances of events 
to the major axis (D). 

All models show a strong excess of far-disk events, with A increasing from 5.3 to 7 as the models go from 
tangential to radial anisotropy. It is also evident from D in Figure 1 that far-disk events lie systematically farther from 
the major axis than near-disk events, providing a second signature of asymmetry. However D is nearly constant across 
the range of a so this spatial signature does not provide a probe of the degree of velocity anisotropy. 

There is also an asymmetry in the Einstein crossing times of M3 1 MACHOs. If MACHOs have radially distended 
orbits, then their motion tends towards being parallel to the line of sight for the near disk but orthogonal to it for the 
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far disk. Events on the near side therefore last longer. In other words, the ratio of the typical time-scales of far-disk 
to near-disk events decreases with increasing radial anisotropy. Consequently, the ratio of the numbers of events in 
the far-disk to the near-disk is enhanced in radially anisotropic models compared to isotropic model. This can be seen 
from the fact that the spatially-averaged microlensing optical depth, (r) oc (?E)(rp), is independent of the velocity 
distribution. From Figure 1, the expected enhancement in the number asymmetry A is ~30%. 

For the isotropic model (a = 0), the time-scales are typically shorter in the far disk than in the near disk. The 
reason for this is that the typical separation between lens and source is larger for near-disk MACHOs. For far-disk 
events, the typical lens-source separation is biased towards the location where the density peaks along the line of sight, 
which is at a distance |j|tanj in front of the sources, where y is the projected distance along the minor axis. The 
situation is different for near-disk MACHOs, where the line of sight density is always a monotonically decreasing 
function of lens-source separation. Here, the typical separation is ~20-30 kpc for y > 5', so {tE)i/ (?E)n is less than 
unity for an isotropic model, as shown in Figure 1 . 

Unfortunately, the Einstein crossing time is not generally measurable for pixel-lensing events. Instead experi- 
ments measure the full- width half-maximum (?i/2) of the lightcurve, which is additionally correlated with the source 
luminosity and background surface brightness distribution. The resulting f|/2 distributions are therefore predicted to 
be broad and any asymmetries in the time-scale distribution are unlikely to be easily observable. 

3. The Effects of Flattening 
3,1, Flattened Models 

Self-consistent solutions of the self-gravitation equations for the density, potential and velocity distributions of 
flattened halo models are rare. Even solving the Jeans equations for the second velocity moments can lead to cumber- 
some results. Baltz, Gyuk & Crotts (2003) have computed pixel-lensing rates explicitly for simple analytic flattened 
halo models. However, if we are merely interested in comparing the effects of flattening on the ratios of quantities in 
the near- and far-disk, then there is a quick alternative to carrying out computations with a fully axisymmetric halo 
model. 

Figure 2 shows lines of sight passing through an elliptical halo with axis ratio q and eccentricity (l-q^Y^^. The 
lines strike the disk of M31 at an angle i. The ratios of the optical depths in the near and far disk are proportional to 
the line segments NP and FP. Also shown is an equivalent spherical halo. By extending N vertically to N' and F to F', 
we can construct lines of sight that pass through the spherical halo and strike the disk at a different angle /'. Similarly, 
the ratios of the optical depths are proportional to the Une segments N'P and F'P. From the elementary properties of 
the elUpse, it follows that the two ratios are the same (NP/FP = N'P / F'P). By straightforward trigonometry, one has 
<7tan; = tan/'. In other words, the asymmetry signal of a flattened halo with axis ratio q is the same as that of an 
equivalent spherical halo, provided the disk is viewed not at angle / but at an angle /'. 

This transformation takes into account the geometric effects of flattening. The first-order changes in the velocity 
distribution can be computed using the tensor virial theorem (Binney & Tremaine 1987; Han & Gould 1995). Strictly 
speaking, the tensor virial theorem apphes globally and relates the components of the total kinetic energy tensor T to 
the components of the total potential energy tensor W. If we assume that the virial theorem holds at each spot, then it 
follows that 

^«Zf£ = !?««l. (5) 

w^z q^' 

Here, we have assumed that the figure is oblate spheroidal, with the the short axis being in the z direction and the 
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plane being equatorial. Although this is only an approximate relation, valid for small flattening, it shows that the 
first-order changes in the velocity distribution are also accounted for by the transformation. 

This means that a quick way to study the asymmetry properties of flattened haloes is to take isotropic spherical 
models and vary the inclination angle of the M3 1 disk. We stress that the transformation does not allow us to calculate 
absolute quantities like the rate, but only the ratios of such quantities in the near and far disk. 



3.2. Results 

Figure 3 shows the variation in the far-to-near-disk ratios for the pixel lensing rates A, the means of vertical dis- 
tances D and the time-scales. The asymmetry signal A for a halo with flattening q is linearly related to the asymmetry 
signal for a spherical halo Aq, that is 

A~l+^(Ao-l). (6) 

The asymmetry signal clearly diminishes with flattening; it is obvious that, in the completely flat limit where the halo 
becomes a razor-thin disk, the asymmetry must vanish. The change in A with flattening is almost entirely caused by the 
change in the ratio of optical depths. Therefore the ratio of the average time-scales is largely unaffected by flattening. 

The distance asymmetry D also decreases with increasing flattening, as is again obvious in the razor-thin Umit. 
The distance asymmetry arises in a spherical model because of two effects. First, lines of sight are longer as we move 
from near to far disk. Second, the line of sight with greatest column density goes through the centre in the near-disk 
but lies away from the centre in the far-disk. So the distribution of distances of events is monotonically decreasing 
in the near disk, but rises to a maximum and then decreases in the far disk. As the flattening increases, all Unes of 
sights become shorter and the density becomes more concentrated towards the centre. The distributions of distances 
in both the near and far disk shrink and the maximum moves towards the centre in the far disk. This latter effect is the 
dominant one, and so the distance asymmetry signals faUs with increasing flattening. 



4. Signal Detection 

4.1. The Number Asymmetry Signal 

In this section, we ask whether one can detect the asymmetry signals in the presence of contaminating events and 
how the velocity anisotropy and the halo flattening affect this detectability. That is, we ask how many candidates are 
needed to detect the asynnmetry at a certain confidence level. 

First, we consider the number asymmetry A between the far and near side. One problem is that none of the 
experiments will be able to obtain pure samples of M31 MACHO events; some microlensing contamination from 
Milky Way MACHOs and M3 1 stars, as well as non-microlensing contamination from variable stars and supernovae 
in background galaxies will be inevitable. Contamination by periodic variables can be minimized by observing over 
a sufficiently long time. Colour information can also be used to eliminate variable stars. Fortunately, the other 
contaminants are equally likely to occur in the near and far regions of the M3 1 disk. Whether or not an asymmetry can 
be detected is therefore a question of the size of the sample, the magnitude of the underlying M3 1 MACHO asymmetry 
and the level of contamination. 

Let us denote the numbers of M3 1 MACHO events in the near and far disk by A'n and A'f respectively, and the 
number of contaminants by A^c. The condition for a detection of the asymmetry signal at the s-a level is that the 
difference between the far and the near counts be greater than s times the Poisson error, which is given by the square 
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root of the total number of the events. Therefore 



y/Nf + Nn + Nc 



>s 



(7) 



where A^j' and A'^ are the total number of candidates on the far and near sides, respectively, including contaminant 
events. In the case where the contaminants are distributed evenly between the far and near disk, A^^ -A^^ = Nf-Na. 
Straightforward manipulation of the above condition leads to a condition on the total number of candidate events, 



where A = Nf/N^ is the underlying asymmetry of the M31 MACHO events, /c = Nc/iN{+Nn) is the contamination 
factor, and A'i = Nf+Nn+Nc the total number of candidate events. Figure 4 shows the expected number of events 
(including the contaminants) required to detect an asymmetry at the 99% confidence level (s = 2.58) for various 
contamination factors. We see immediately that the size of the asymmetry signal is crucial. If A < 2, then even 
with low contamination the number of M3 1 microlensing events needed to give a convincing detection of asymmetry 
exceeds 60. On the other hand, if the asymmetry is ~5, then more than 15 candidates are required. Reassuringly, on 
referring to Figures 1 and 3, we see that A typically Ues between 5 and 7 for spherical models, and only becomes as 
small as 2 for models with q = 0.3. In fact, haloes flatter than this are not likely on dynamical grounds, as they are 
susceptible to bending instabilities (e.g., Merritt & Sellwood 1994). Samples of ~15 candidate events are well within 
reach of the current surveys if MACHOs contribute significantly to the dark matter budget. However, this number is 
only the contribution of the M31 halo, as we have cut out the central parts of the M31 disk. 

If the contribution of MACHOs to the dark matter mass budget is significant in M3 1 , we expect a similar MACHO 
contribution for the Milky Way. For interesting MACHO fractions, the foreground Milky Way MACHOs may well 
provide the dominant contribution to f^. If the typical MACHO mass and halo density contributions are universal, 
then the magnitude of fc will be determined, to first order, by the relative masses of the M31 and Milky Way haloes. 
If the M31 halo is twice as massive as the Milky Way's, then /c ^ 0.4 (Kerins et al. 2001), whilst we should expect /c 
to be closer to unity if the two haloes are equally massive (Evans & Wilkinson 2000). In any case, the total number 
of candidates required to confirm asymmetry scales as (1 + /c)^ from equation (8), so we require four times as many 
candidates to detect the asymmetry when /c = 1 as when = 0. 



Motivated by the variation in mean distances of events between the near and far disk seen in Figures 1 and 3, we 
can go beyond the asymmetry in number counts by considering a distance asymmetry signature. 

To quantify differences between the near-disk and far-disk event positions, we apply the Mann- Whitney rank-sum 
test (Mann & Whitney 1947). This is a non-parametric test for differences between the medians of two samples. The 
Mann- Whitney test exploits the fact that two samples drawn from identical distributions exhibit the property that, if 
one combines them and then ranks the elements by size, the two samples are uniformly intermingled on average within 
the ranked combined sample. For two samples of large enough size and rib, the sum of the rank numbers for sample 
a is normaUy distributed about a mean Ha = naina+rib) /2 with a variance = naftbina+rib+V) / \2. The nuU hypothesis 
of similarity (or, in our case, symmetry) is therefore straightforward to quantify. 

In the case of M3 1 pixel lensing, we can apply the Mann- Whitney test to the distribution of y, the projected 
distance from the M3 1 major axis, for candidates in the near- and far-disk sub-samples. We can also combine the 




(8) 



4.2. The Distance Asymmetry Signal 
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Mann-Whitney test with the number asymmetry test, since the Mann- Whitney statistic probes the spatial distribution 
of events, not the relative sizes of the samples. If the far-disk event positions are designated as sample a and the sum 
of their rank numbers within the combined near- and far-disk sample is 6a, then we can define ^mw = (^a~A*a)/ca- 
Taking to be the significance of the number asymmetry statistic, as defined by the left-hand side of equation (7), 
then the overall significance of the combined sample iss = (i^^^+i^)^^^- A value ofs = 2.58 would indicate a near-far 
asymmetry favoured at the 99% confidence level. 

We have performed Monte Carlo simulations to determine the number of candidates required to secure a 99% 
confidence detection of asymmetry. The simulations test a range of anisotropic and flattened halo models and a range 
of contaminations. For a given model, event position realisations are generated using the theoretical pixel-lensing rate 
of equation (3), weighted by source number density. For each new event realisation, the overall significance, s, of 
the cumulative sample is computed from both the Mann- Whitney statistic (^mw) and the number asyrmnetry statistic 
(sn). When there are fewer than five events in either the far- or near-disk sub-samples, only the number asymmetry 
statistic is used because the distribution of rank sums for small data-sets can deviate strongly from Gaussianity. We 
assume the contaminating populations are symmetrically distributed about the M3 1 major axis, though we adopt the 
most difficult case when they have a comparable spatial dispersion to the M3 1 MACHO events. The extent to which 
this is true depends upon whether Milky Way MACHOs, on the one hand, or variable stars and stellar microlenses, 
on the other hand, provide the dominant contribution to N^- If the latter population dominates, then contaminants are 
likely to be more spatially concentrated than M3 1 MACHOs, making the MACHO asymmetry easier to measure for 
a given /c. So for each pixel-lensing realisation, a uniform random number in the interval [0,1] is chosen. If it is less 
than /c/[2(l +/c)], then the position of the event is flipped about the M31 major axis. This means that, on average, 
a fraction /c/(l + /c) of the total sample is symmetrically distributed, as required. A full trial, is terminated either 
when the cumulative sample of A^t,; candidates provides at least a 99% confidence detection of asymmetry, or when 
the program estimates that the required sample is likely to exceed 500 events. This whole process is repeated 1000 
times for each model and the median value of A't,, (excluding trials which are prematurely terminated) is adopted as 
the estimate for Nt. 

The upper panel of Figure 5 shows the median number of candidate events (including the contaminants) required 
to detect an asymmetry with 99% confidence, as a function of the halo velocity anisotropy parameter a. In the absence 
of contamination, a sample size between 1 1 to 14 events is typically sufficient to detect asymmetry. The raggedness of 
this line is partly due to Monte Carlo noise but also partly due to the fact that A't necessarily takes only a few discrete 
values in the limit of small data-sets. From equation (7), the smallest sample needed give rise to a 99% confidence 
detection of asymmetry is 7 events, all of which must be located in the far disk (assuming the asymmetry is caused 
by M31 MACHOs). When one candidate lies in the near disk we need A'i = 11, and when two candidates lie in the 
near disk we require A^t = 14. The leaps in the thick solid line between 1 1 and 14 events reflect this discrete behaviour, 
though the oscillation back and forth for -2 < ln(l + a) < -1 is due to Monte-Carlo noise. In the worst case considered 
in Figure 5, where contaminants outnumber M31 MACHOs by 3:1, a median sample of ~ 120 events is required to 
detect asymmetry. For a = and /c = 3, we find that the expectation value of A'i is around 190. Using s = 2.58, /c = 3 
and an asymmetry A = 5.5 when a = (Figure 1), equation (8) indicates that we should require A^t = 220 if we use 
number-count information alone. Therefore the addition of distance information aUows around a 15% reduction in the 
required number of candidates in this case. 

Overall, A^t does not appear to be particularly sensitive to a. This is to be expected when contamination levels are 
high because the 30% contrast in A between the radiaUy- and circularly-anisotropic models (see Figure 1) is strongly 
diluted. 

The lower panel of Figure 5 shows the situation for flattened halo models. The larger range in A for the flattened 
models (1 .3 < A < 5 .5) means that the median Nt shows a greater sensitivity to flattening than to velocity anisotropy. 
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In the absence of contaminants, A^t ^ 25 is typically needed to confirm asymmetry for models rounder than q = 0.3. 
A sample of 100 candidates would permit asymmetry to be detected even if /c — 1. Adding distance information is 
particularly effective at reducing the required size of candidate samples for highly flattened models. As an extreme 
example, when ^ = 0.1 and /c = 0, the expectation of Nt is around 80 if distance information is used along with number 
counts. From equation (8), 750 candidates are required if only number counts are used. This shows the value of the 
Mann- Whitney statistic. 

It is of course much more difficult to measure an asymmetry than to detect one. Suppose an experiment has 
gathered M ^ 100 events (including contaminants). With reference to Figure 4, if there has been no detection of 
asymmetry at the 99% confidence level, then - as A > 3 for all models we have considered with q >0.5 - we can 
infer that the signal has been overwhelmed by contaminants (/c > 1). The degree of contamination must be greater 
than a critical value which is given by computing the curve which passes through the point (A = 3,A^t). So, a null 
signal can be used to give a constraint on the contamination. On the other hand, if there has been a detection, then 
the measured signal A' ~ iNi+jNc)/iND+^Nc) is merely a lower limit to the true signal. It should be possible to 
estimate the contamination fraction statistically even if we do not know the individual contaminating events, so the 
true asymmetry signal can be matched to models using standard Bayesian hkelihood methods (Kerins et al 2001). 
Current surveys should therefore be able to discriminate between halo models with different degrees of flattening 
and may, if the spatial distribution of contaminant populations is well characterized, be able to distinguish between 
radially-anisotropic halo models and tangentiaUy-anisotropic or isotropic models. 



4.3. Confusion from M32 and its Stream 

Additional confusion of the asymmetry signal may come from microlensing by stars which belong to streams or 
tidal debris from disrupted sateUite galaxies cannibalized by M31 and/or by stars belonging to the intervening dwarf 
elUptical M32. 

Ibata et al. (2001) have traced out a giant stream in the M31 stellar halo in red giant branch star counts. The 
stream is ~1° wide in projection. It seemingly originates from the satellite galaxy M32 and possibly also incorporates 
NGC 205. The average surface brightness of the stream is ~30 mag arcsec"^ in the V-band. It has been suggested 
by Ferguson et al. (2002) that this may confuse the detection of the near-far disk asynmietry in the microlensing 
experiments. For sources in M31 and lensing populations at roughly the same distance from the sources, the optical 
depth of the stream is 

T ~ 4.3 X 10-1" ( ( -^^^) 10'2-"-4^ (9) 



v20kpcy VM©/L0/ 

where fi is the surface brightness in magnitude per arcsec^, d the separation of the stream from M3 1 disk along the 
line of sight and M/L is the MACHO mass-to- light ratio of the stream. This ratio must exceed 100 for the optical 
depth of the stream to be comparable to the other lensing populations. Accordingly, the presence of such streams is 
not likely to be problematic for asymmetry detection. 

More problematic may be the intervening dwarf elUptical galaxy M32. The POINT-AGAPE collaboration has 
already found a candidate which lies ~3' in projection away from the centre of M32 and argued that the lens most 
probably lies in M32 itself (Paulin-Henrikkson et al. 2002). The optical depth of M32 is estimated to be (PauUn- 
Henrikkson et al. 2003) 

r^lA><lO-^(^](J^] (10) 



V20kpcy V3M0/L0/ 

where d is the separation of M32 from M3 1 along the fine of sight and M/L is the stellar mass-to-light ratio of M32 
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This is comparable to the signal expected from the baryonic hiilo (assuming a 20% MACHO fraction). Although 
microlensing associated with M32 is a potentially significant contaminent the affected region can be readily excluded 
from the statistics using an appropriate mask. 



5. Conclusions 

Pixel-lensing experiments targeting M31 are hoping to exploit the favourably high disk inclination in order to 
detect an asymmetry in the spatial distribution of microlensing events. If such a signjil is found, it will provide 
powerful evidence for the existence of MACHOs. 

For a spherically symmetric and isotropic halo, the numbers of M3 1 MACHOs in the far-disk outnumber those 
in the near-disk by more than 5 to 1. This asymmetry is increased by about 30% if M31 MACHOs occupy radial 
orbits rather than tangential ones. The signal is diminished if M31 MACHOs lie in a flattened halo. However, even 
for haloes as flat as q = 0.3, the numbers of M31 MACHOs in the far side will outnumber those in the near side by a 

factor of ^2. 

The key to detecting asymmetry is to isolate microlensing events solely due to M3 1 MACHOs. There is likely 
to be significant contamination from other microlensing populations, as weU as from variable stars and supemovae 
mistaken for microlensing, which wiU dilute the observed signal. The combination of number-count and distance 
information permits asymmetry to be detected for a wide range of halo models, even in the presence of significant 
levels of contamination. For models with high levels of asymmetry, such as spherical haloes or haloes with a high 
degree of radial velocity anisotropy, number count information alone provides a sensitive diagnostic. The addition 
of distance information aUows ^15% reduction in the size of samples needed to confirm asymmetry. For models 
with low levels of asymmetry, such as strongly flattened haloes, distance information can reduce the required size of 
candidate samples by a factor 2 or more. 

A sample of 50 events is typically sufficient to detect asymmetry in the M3 1 MACHO distribution within spherical 
haloes, even if only half the sample is due to M31 MACHO events. For flattened halo models, a sample of 50 
candidates would likely allow asymmetry to be seen, provided that the halo axis ratio ^/ > 0.3 and the contaminants do 
not contribute more than a third of the sample. The term "contaminants" covers Milky Way MACHOs, M3 1 disk stellar 
lensing events and variable star populations, all of which are assumed to be symmetrically distributed with respect to 
the major axis. For comparison, Paulin-Henriksson et al. (2003) have already found 362 lightcurves compatible with 
microlensing from the first two years of the POINT- AGAPE survey, though the contamination factor may still be 
very large. Samples of 50 events with modest contamination are easily achievable with the current generation of 
pixel-lensing surveys. 

NWE and JA thank the Royal Society and the Leverhulme Foundation respectively for financial support. 
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Fig. 1. — The spatially-averaged pixel-lensing rate (Fp) {solid line) as a function of anisotropy parameter a and 
normalised to the value for the isotropic model (Fo) = 7.6 x 10"^ stars"^ year"^ Also shown is the M31 MACHO 
number asymmetry A (dashed line), the ratio of near-disk to far-disk average durations {t-Eji/ (?E)n {dot-dashed line), 
and the ratio of projected distances to the major axis D {dotted line) 



-12- 




Fig. 2. — This shows Unes of sight through a spherical and an elliptical halo. From the properties of an ellipse, we 
know that the ratios NP:FP and N'P:F'P are equal. This enables us to relate the asymmetry signal of a disk viewed 
through a flattened halo at incUnation i to the same disk viewed through a spherical halo at incUnation i'. 
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Fig. 3 . — The far-to-near-disk ratio for the pixel lensing rate A {dashed line), the mean vertical distance D {dotted line) 
and the time-scales {dot-dashed line), shown as a function of flattening. This diagram is drawn using the transformation 
introduced in Section 3.1. 




A = N,/N^ 

Fig. 4. — The number of candidate events (M3 1 MACHOs and contaminants) required to confirm asymmetry with 
99% confidence, based upon number-count information alone, plotted as a function of the underlying M3 1 MACHO 
number asymmetry A. The different lines correspond to contamination factors, /c = {thick solid line), 0.3 {thin solid 
line), 0.5 {dashed line), 1 {dot-dashed line) and 3 {dotted line). 
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Fig. 5. — The median number of candidates (M31 MACHOs and contaminants) required for a 99% confidence de- 
tection of asymmetry. Line coding is the same as in Figure 4. Upper panel shows Ni as a function of halo velocity 
anisotropy for a spherical halo. Lower panel shows Nt for flattened halo models, using the equivalence between the 
halo flattening parameter q and disk inclination ; discussed in Section 3.1. 



